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^ ' Abstract 

Nuclear-polar izability corrections that go beyond unretarded-dipole approx- 
imation are calculated analytically for hydrogenic (atomic) S-states. These re- 
tardation corrections are evaluated numerically for deuterium and contribute 
-0.68 kHz, for a total polarization correction of 18.58(7) kHz. Our results are in 
agreement with one previous numerical calculation, and the retardation correc- 
tions completely account for the difference between two previous calculations. 
The uncertainty in the deuterium polarizability correction is substantially re- 
duced. At the level of 0.01 kHz for deuterium, only three primary nuclear 
observables contribute: the electric polarizability, aE, the paramagnetic sus- 
ceptibility, Pm, and the third Zemach moment, {r^)(2)- Cartesian multipole 
decomposition of the virtual Compton amplitude and its concomitant gauge 
sum rules are used in the analysis. 
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Introduction 



The remarkable experiments presently being performed in Garcliing|]I| and ParisjQ 
on the spectroscopy of hydrogen isotopes have astonishing precision. The Rydberg 
currently has an uncertainty of 9 parts per trillion, while the isotope shift between 
deuterium and hydrogen IS - 2S transitions has a reported uncertainty of 3 parts per 
billion, and this is expected to be lowered soon by an order of magnitude0. 

The isotope-shift measurements afford a unique opportunity for nuclear physics. 
The traditional technique for determining nuclear sizes is to scatter relativistic elec- 
trons from nuclei, determine the charge form factor, and extrapolate this to small 
momentum transfers, thus determining the mean-square charge radius, (r^). It is 
extremely difficult to perform the latter measurements with an absolute accuracy of 
one percent or less, and this sets limits on the accuracy of the charge radius. The 
currently accepted value of the charge radius of the protonQ, (r^)^^ = 0.862(12) fm, 
corresponds to an uncertainty in (r^) of nearly 3 percent, and the recently determined 
deuteron radius0, (r^)y^ = 2.128(11) fm, has an uncertainty in (r^) of 1 percent. For 
the 1S-2S d-p isotope shift [|ll the nuclear-size correction contributes approximately 
-5000 kHz (roughly the same as the QED corrections) out of a total of 670 GHz. The 
reported]^ uncertainty of 2 kHz corresponds to a precision of better than one part 
per thousand in (r^). 

In addition to static size corrections, the electron polarizes the nucleus and pro- 
duces nuclear-polarizability corrections. In order to use the isotope shift as a precise 
gauge of nuclear size differences, it is necessary to compute these polarizability cor- 
rections as accurately as possible, and that is the goal of this work. 



There have been several calculations of these corrections for deuterium |Q |^, ||, H, 



10, 11|. The bulk of the effect (^19 kHz in toto) is caused by the Coulomb interaction 
distorting the nucleus (~17 kHz) with a smaller (;^2 kHz) contribution from the 
virtual transverse photons. In leading order (unretarded-dipole approximation) the 
electric polarizability, as, dominates the process and accounts for 19.26(6) kHz in 
nonrelativistic approximation for the deuteron |Tl| . This numerical result summarized 
calculations for a group of "second-generation" potentials, which fit the nucleon- 
nucleon scattering data well enough to be considered alternative phase-shift analyses 
of that data. 



There exists a single calculation 0] , using first-generation nuclear potentials, that 
goes beyond the unretarded-dipole approximation and includes retardation, higher 
multipoles, the effect of the finite sizes of the nucleons, seagulls, and even meson- 
exchange currents. Results for this calculation are smaller by > 0.5 kHz than for 
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Figure 1: Nuclear polarization corrections with direct (a), crossed (b), and seagull 
(c) contributions are illustrated. Single lines represent an electron, double lines a 
nucleus, and shaded double lines depict an excited nucleus, with the seagull vertex 
maintaining gauge invariance. 

those using the unretarded-dipole approximation. That calculation was performed 
by constructing nuclear charge and current (transition) densities and performing a 
difficult double integral over the momentum and energy transferred across each pho- 
ton line (see Figure 1). Our goal is to reduce that calculation to an analytic series in 
various size- dependent nuclear observables, and keep only those that are expected to 
contribute at a level of greater than 0.01 kHz. The resulting expression is fairly simple 
and depends only on three primary nuclear observables: the electric polarizability, ue-, 
the paramagnetic susceptibility, (3m, and the third Zemach moment [jl2|, |13|, (r^)(2) , 
of the charge distribution. Meson-exchange currents |ll4| play a small role that is eas- 
ily incorporated in the calculation. Our final result is in excellent agreement with the 
difficult but comprehensive calculation of Ref. 0. We will produce a final estimate 
for the complete polarizability correction of 18.58(7) kHz, based on second-generation 
potentials. It will be not be easy to improve this result significantly, because it will 
be difficult to increase the precision of the nuclear observables on which the result 
depends. 



Higher Polarizabilities 



The integrals over momentum transfer in the loops that define the generalized 
polarizability correction are difficult, rather complicated, and extremely tedious. For 
all these reasons, we have relegated them to Appendix A. The constraints of gauge 
invariance are crucial to impose (the results are infrared divergent otherwise), but are 
also tedious to develop, although they have been known for decades[|l5| . Consequently, 
a brief presentation of the necessary relations has been relegated to Appendix B. Only 
those parts of the calculation that we will treat numerically are given directly below. 
To the order that we work, gauge invariance has been properly implemented. 
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We first define the electric polarizability||T5|, [T6|, a^;, in terms of the electric-dipole 
operator, D, 

«E = ^2^-7^ 7^, (1) 



and the logarithmic mean-excitation energy [pT] , p!7| , ii^, by 

log(^/me) ai, = ^ E ^^^^ log [(i?^ - i?o)/me] . (2) 

Similarly, we define the paramagnetic susceptibilityll^, |l6l, [3m , in terms of the 
magnetic- (dipole) moment operator, /x, (see Eq. (B8)) 

_2a^ KiV|M|0)p 
Pm - — 2^ -7^ 7^ , l<3j 

together with its mean-excitation energy, 

log(E/me) /?M = ^ E 'If' log [(^^ - ^o)/m,] . (4) 

A close relative of a^; is 

2a |(iV|D|0)p 

together with its logarithmic mean-excitation energy analogously defined. Although 
we have used (and will write below) E for each of the different mean-excitation 
energies for notational simplicity, they are distinct (although similar) numbers. They 
will always be grouped with the operators that define them. In the formulae above, 
a is the fine-structure constant, me is the electron mass, |A^) is the Nth nuclear state 
(the Nth eigenstate of Hq) with energy, Ej^, and iV = labels the ground state. 

The inelastic charge density (squared) can be rewritten in cases where there are 
no energy factors] 



E(0|pt(x)|Ar)(Ar|p(y)|0) = (0|pt(x)p(y)|0) - po(x)po(l/) , (6) 

where 

Po(x) = (0|p(x)|0), (7) 
p(x) is the nuclear charge (-density) operator and / dPx po{x) = Z, the nuclear charge. 
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With these definitions, we can rewrite Eq. (A30) in the form 



AEn = -4ame|(/)„(0)p 



5aE 



19, 15 



283, 



[log(2E/me) + — 1 + — Aas[log(2E/me) - 
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+ 
where 



[log(2E/me 



1 7ia 
6^ ^ T2 



d'x / A|x-y|' (0|pt(x)p(y)|0)-po(x)po(z/) 



A5 = + log (Ez) - ^) I d^a; | d'y (0| Afi(x, y)|0) 



and 



Ai?(x, y) = S^^(x, y)[x . y - y^] + BHi^, y)[-2yy 



(9) 



(10) 



We have dropped a large variety of small terms in AB (indicated by dots and with 
relative coefficients ~1/10) that arise from the seagull and current terms (converted 
using gauge sum rules; see Appendix A). Gauge sum rules and approximating E^ — Eq 
by E and 2; = |x — y| by ^ (i.e., constants) were used to obtain Eqs. (9) and (10). 
The small correction AB arises from the Coulomb interaction; we will also use the 
alternative expression in Eq. (A31). Note that the term ~ |x — y|^ contributes to all 
electric multipoles, unlike the others, which are either electric- or magnetic-dipole in 
nature. Only ground-state properties are needed to construct the former term. 

In order to proceed further, we need to specify our nuclear model. The charge 
operator p(x) is given by 



P(x) = ^ei(x-Xi) 

i=l 



where 

Ci = ep{i) ep(x - x^) + e„(i) e„(x - Xj) (12) 

counts protons (with e^) and neutrons (with e„) and multiplies each species by its 
intrinsic charge distribution. The form (11) is nothing more than the usual folding of 
ei{x) with (5^(x — Xj). Forming / ci^a:xp(x) = I]iep(i)xj demonstrates that finite size 
does not modify the electric-dipole operator. 

The current operator needed to construct n consists of three distinct parts: the 
spin-magnetization current, the orbital current, and meson-exchange currents (MEC). 
We ignore nucleon finite size, which does not contribute to this order, and find[|T^] 



i=i \ 



2M 



f^MEC ■■ 



(13) 
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where the spin-magnetization current is determined by 



fl{i) = fipCpii) + flnCnii) ■ (14) 

Note that the isoscalar and isovector nucleon magnetic moments are very different in 
size: fig = f^p + fJ'n = 0.8798 . . . and /i,; = /ip — /U„ = 4.7059 . . .. The large value of the 
isovector nucleon magnetic moment will play a determinative role. We eschew writing 
out our model for the two-body pion-exchange currents(i.e., MEC), which is discussed 
in Ref. []TU|. This model has had its pion-nucleon form factor adjusted to reproduce 
the experimental thermal n — p radiative capture rate pO. As the contribution of 



is relatively small and the MEC a small part of this, the overall MEC contribution is 
nearly negligible, but has been included for completeness. 

Our final ingredient is the Compton seagull operator. This operator is comprised 
of several component spT|: impulse approximation, plus meson-exchange currents [ p^] , 
plus .... We expect the meson-exchange currents to be possibly comparable to the 
impulse approximation, based on sum-rule studies|]23|. We will work only with the 
impulse- approximation component, which has the form[|l^ 

A 

5*^(x, y) = ^ e,(x - X,) e,(y - x,) . (15) 

2=1 

The pion-exchange component of the deuteron's diamagnetic susceptibility (see Eq. 
(A29) and Refs. [0, ^) has been shown to be tiny. 



Numerical Calculations 



Podolsky's method|2^ is very convenient for calculating as and Pm- Any gen- 
eralized polarizability of the type displayed in Eqs. (1) and (3) can be calculated as 
follows. Equation (1) is fully equivalent to 

aE = 2a{0\D,\A^,), (16) 

where 

(Ho-Eo) |A^,) = D,\0) (17) 

is solved subject to finite boundary conditions. One must be careful to exclude the 
ground state from the sum over N in Eq. (3) for /Sm- This necessitates a projection 
orthogonal to the ground state on the right-hand side of Eq. (17) (with replaced 
by fj,z in both Eqs. (16) and (17)). 



For the deuteron, one impulse- approximation calculation of (3m exists[^ with a 
value of 0.065 fm^. This is dominated by ^5*0 intermediate states. Indeed, an upper 
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limit for all triplet intermediate states is (3j^j < 0.0003 fm^, obtained[^ by replacing 
the energy denominator in Eq. (3) by its smallest possible value {Ed, the deuteron 
binding energy), and then using closure and completing the algebra. Moreover, the 
^D2-state contribution is tiny, and the ^5*0 intermediate states dominate completely. 

The logarithmic mean-excitation energies are calculated using a trick|2^. We 
define a quantity closely related to 

MO^IeJ™^. (IS) 

where / is an energy-scaling factor ~ 3-5-5'^ inserted for convenience. An integral 
over ^ produces a logarithm, and one finds a convenient numerical algorithm for E in 

as(0) log (2E/me) = /' - «£;(0) + ^^(l/O] - «i^(0) logK/2/) , (19) 

Jo 4 

where E is independent of /, and Eq. (19) is fully equivalent to Eq. (2). 



The electric polarizability was calculated and thoroughly discussed in Ref. ITT 
One found there 

aE = 0.6328(17) fm^ 
log(2^/me) = 2.9620(5) (20) 
z/p„^ = 19.26(6) kHz = (16.98) + (2.28) kHz , 

where the latter is broken down into Coulomb, transverse, and total contributions. 
This calculation did not incorporate relativistic corrections to the deuteron. 

The analogous calculations for Pm are listed in Tables I and II for a variety of 



first- and second-generation potentials [26, U7\, Eq, B3, BO, 31, 32, 33, p3, pH]. The 



results of Table I are for the impulse-approximation magnetic moment (no MEC), 
while those of Table II incorporate MEC, as well. The latter increases the former by 
approximately 15%, which is quite typical for isovector transitions. We average the 
various second-generation results and estimate 

= 0.0777(3) fm^ 
log(2^/me) = 2.498(2) (21) 
i/p^^i = -0.307(2) (6) kHz. 

These uncertainties do not include uncertainties in the MEC, which are possibly 1- 
2% of the total result (this is subjective); the latter is refiected in the second error 
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Table 1: Impulse-approximation deuteron magnetic susceptibilities, /5m, in units 
of fm^, logarithmic mean-excitation-energy ratios, log(2^/me), and corresponding 
deuteron 1S-2S polarization-energy shifts, z/poi, in kHz. The RSC potential labelled 
[-23] has had its ^Sq part modified to produce the correct n-p scattering length. 



Potential Model 




\og{2E/m,) 


z/poi(kHz) 


Second-Generation Potentials 


Argonne Vis 


0.0678 


2.4724 


-0.265 


Nijmegen (loc-rel) 


0.0677 


2.4726 


-0.264 


Nijmegen (loc-nr) 


0.0677 


2.4732 


-0.264 


Nijmegen (nl-rel) 


0.0677 


2.4726 


-0.264 


Nijmegen (nl-nr) 


0.0676 


2.4732 


-0.264 


Nijmegen (full-rel) 


0.0675 


2.4728 


-0.264 


Reid Soft Core (93) 


0.0674 


2.4744 


-0.264 


First-Generation Potentials 


Bonn (CS) 


0.0682 


2.4738 


-0.267 


Argonne V14 


0.0674 


2.4733 


-0.263 


Reid Soft Core (68) [-23] 


0.0669 


2.4748 


-0.261 


Nijmegen (78) 


0.0663 


2.4947 


-0.261 


Super Soft Core (C) 


0.0659 


2.4982 


-0.260 


de Tourreil-Rouben-Sprung 


0.0656 


2.4969 


-0.259 


Paris 


0.0653 


2.5008 


-0.258 


Reid Soft Core (68) 


0.0647 


2.5031 


-0.256 



(6) in the last Eq. (21). This polarization contribution is nonnegligible only because 
fly = 4.7; a more "normal" size ~1 would reduce the contribution by a factor of ~25. 
Our z/poi is in reasonable agreement with the zero-range result of Ref. [^. 

The very small corrections, Actg, can be accurately estimated in zero-range ap- 
proximation (which we used as a tool in Ref. [ll| ) . We first calculate (^ ) using Eq. 
(18) and find 

where the asymptotic (reduced) s-state wave function of the deuteron has the form: 
As e~'^^. Moreover, R,"^ = fi;^(l + ^/), and /i is the n-p reduced mass. Performing two 
derivatives leads to 

ml Aa% = ^-^f^ = 0.0098 fm^ , (23) 



24E 



and is accurate to better than 1/2%. The following integral (similar to Eq. (19)) 
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Table 2: Full deuteron magnetic susceptibilities, (3m, in units of fm^, logarith- 
mic mean-excitation-energy ratios, log(2^/me), and corresponding deuteron 1S-2S 
polarization-energy shifts, z/poi, in kHz. The RSC potential labelled [-23] has had its 
^Sq part modified to produce the correct n-p scattering length. 



Potential Model 




lGg(2E/m,) 


;/poi(kHz) 


Second-Generation Potentials 


Nijmegen (fuU-rel) 


0.0780 


2.5003 


-0.308 


Nijmegen (nl-rel) 


0.0778 


2.4981 


-0.307 


Nijmegen (nl-nr) 


0.0777 


2.4987 


-0.307 


Nijmegen (loc-rel) 


0.0775 


2.4972 


-0.306 


Nijmegen (loc-nr) 


0.0774 


2.4978 


-0.306 


Reid Soft Core (93) 


0.0775 


2.5005 


-0.306 


Argonne Vi8 


0.0774 


2.4963 


-0.305 


First-Generation Potentials 


Argonne V14 


0.0774 


2.4996 


-0.306 


Reid Soft Core (68) [-23] 


0.0769 


2.5002 


-0.304 


Bonn (CS) 


0.0766 


2.4935 


-0.302 


Nijmegen (78) 


0.0751 


2.5172 


-0.299 


de Tourreil-Rouben-Sprung 


0.0748 


2.5221 


-0.298 


Super Soft Core (C) 


0.0747 


2.5218 


-0.298 


Paris 


0.0743 


2.5253 


-0.297 


Reid Soft Core (68) 


0.0742 


2.5295 


-0.297 



produces the logarithmic mean-excitation energy 

1 



log {2E/me) 



I' f [«l(0 - «^(o) - e<(0) - |-<(0)] 

- log(me/2/) 



+ /' d^[^a%{l/0 - e«^(0) - <(0)] 
Jo 

= 2.648. 

These results produce a total correction from AaE' 

Aa 



(24) 



(25) 



i/^J^ = 0.106kHz = (0.060) + (0.046) kHz, 
which has been broken down into Coulomb and transverse parts, respectively. 

The remaining large quantity in Eq. (8) is the retardation correction proportional 
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to pin(x, y). Using Eq. (11), we expand that result and find 

p^(x)p(y) = ei(x-Xi)ej (y-xj)+^ ep(z)ep(x-Xi)ep(y-Xi)+e„(ii)e„(x-x,;)e„(y-Xi 

(26) 

Note that if the neutron's charge distribution is set to zero, the first term vanishes 
for the deuteron (one nucleon must be a neutron), and only the second term survives. 
Shifting the x and y integrals each by Xj in that case produces {r'^Y^2)^ while the 
Po{x)po{y) term generates (r^)l-,, and a retardation correction 



7ra 
I2 



where 



-{r%) + {r%)\^J^Mr')i2), (27) 
(r")(2)= fd\r-p^,^{r), (28) 



and 

P(2) 



(r) = J dhpi\z - r|) p{z) =p0p (29) 

is the convoluted (Zemach^, |l3l) density. In Eq. (27), the first of the third moments 
is calculated with respect to the total deuteron (including the finite size of the proton) 
convoluted charge density and the second with the proton's convoluted density, pp. In 
what follows below we will specialize to the deuteron, and denote by pd the deuteron's 
ground-state charge density (called po before). 

For completeness, we include the neutron contributions as well. There will be 
a two-body correlation term (first term in Eq. (26)) involving p„ (8> Pp Pd, where 
Pd is slightly modified to account for the vector r specifying a correlation, while r/2 
determines the charge density, pd- In addition to this term, the folded proton density 
in Eq. (27) is replaced by {pp® Pp + Pn® Pn) and the deuteron charge density is defined 
by (Pp + Pn) ® where p° is determined by the deuteron wave function alone. 

We use a simplified model of the neutron and proton form factors. The proton 
form factor is taken to have a dipole form with the correct radius^] (0.862 fm). The 
neutron form factor is that dipole times g^, adjusted overall to match the experimental 
charge radius of -0.338 fm|36|. 



Calculations of the various moments are performed by first calculating the deuteron 
density and then generating a spline fit of it. Then folding is performed and that con- 
voluted density is similarly fit. Moments are calculated ultimately with respect to 
the final fitted density. 

The results for various models (including the effect of neutrons) are listed in Table 
3. The neutrons lower the result by approximately 1%. The second-generation results 
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Table 3: Zemach-moment contribution A(r^)(2), to the Coulomb-induced retarda- 
tion correction in units of fm^, and corresponding deuteron 1S-2S polarization-energy 
shifts, z/poi, in kHz. 



Potential Model 


A(r3)(2)(fm=^) i/po,(kHz) 


Sect)ncl-Generation Potentials 


Argonne Vi8 
Reid Soft Core (93) 
Nijmegen (loc-nr) 
Nijmegen (loc-rel) 
Nijmegen (nl-rel) 
Nijmegen (nl-nr) 
Nijmegen (fuU-rel) 


-37.45 -0.485 
-37.44 -0.484 
-37.38 -0.484 
-37.34 -0.483 
-37.36 -0.483 
-37.31 -0.483 
-37.20 -0.481 


First- Generation Potentials 


Super Soft Core (C) 
Nijmegen (78) 
Argonne V14 

de Tourreil-Rouben-Sprung 
Paris 

Bonn (CS) 

Reid Soft Core (68) 


-38.45 -0.498 
-38.27 -0.495 
-38.01 -0.492 
-37.70 -0.488 

-37.55 -0.486 
-37.44 -0.484 
-36.87 -0.477 



for this retardation (and higher-multipole) correction can be summarized by 

A(r^)(2) = -37.32(12) fm^ , (30) 

and 

i/^^t = -0.483(2) kHz. (31) 

This process is the only one that we will consider where higher multipoles, retarded 
El, and nucleon finite size contribute. 

The final task will be to estimate the size of AB. The "natural" size of terms 
with numerical coefficients ~1 is AB ~ while the coefficient 4mea^|0„(O)p 

has a value 0.05 kHz/fm^. Since -^^^ ~ 0.8 fm^, the natural size is 0.04 kHz. We 
use > 10 MeV, z > 4 fm to estimate logarithms and Eqs. (15) and (B12) for the 
seagull operator. The seagull contribution has a size ~ ^ (0.04 kHz) ~ 0.005 kHz. 
The higher-order current terms are of similar size, but largely cancel, leaving a very 

tiny residue. The higher-order charge terms are dominated by quadrupole excitations 
and have a nominal size ~ | (0.04 kHz) ~ 0.04 kHz, which is almost as large as the 



11 



Table 4: Contributions in kHz to the deuteron-polarizability frequency shift for the 
1S-2S transition together with their respective origins, separated into Coulomb and 
transverse, and electric dipole, magnetic dipole, and higher-mult ipole and retardation 
terms. 



Origin 


dE AaE pM A(r^)(2) AB 


Total 


Coulomb 
Transverse 


16.98 0.06 - -0.48 < 0.01 
2.28 0.05 -0.31 - < 0.01 


16.56 
2.02 


Total 


19.26 0.11 -0.31 -0.48 <0.01 


18.58 



AaE Coulomb term. It can be shown, however, that this size is an artifact, caused 
by neglecting some recoil terms. We must perform the estimate more carefully. 

If one uses instead the representation in Eq. (A31) for the z'^ charge terms, and 
evaluates the double commutators, one finds that the kinetic-energy part of Hq van- 
ishes in the point-nucleon limit, and otherwise has a rough size ~ 0.006 kHz. One 
can also evaluate the potential part of the commutator and find ~ 0.003 kHz. These 
corrections are not only small, but the comparable sizes of kinetic and potential terms 
are in accordance with expectations. 



Results and Conclusions 



Our final results are tabulated in Table 4, with breakdowns according to their 
origin. The total for the 1S-2S transition in deuterium is 



poi = 18.58(7) kHz . (32) 



This is 0.68 kHz less than the leading-order {oe) result, and is consistent with the 
differences between the results of Refs. and p. The complete numerical results 
of Ref. 0] for four first-generation potentials (Paris|^, y4Vi4[^], Nijmegen p^ ,and 
Bonn CS|30|) are in agreement with our results within 0.02 kHz for Coulomb and 



transverse parts, which must be regarded as virtually perfect agreement. We note 
that substantially improving the uncertainty in Eq. (32) will be difficult, because it 
would entail substantial improvement s JTT|, P7| in the nuclear parameter, Ag. 



A wide variety of physical mechanisms contribute to the final result. Unretarded 
El photons (both longitudinal and transverse) generate the electric polarizability. 
The paramagnetic susceptibility generates a nonnegligible term only because the nu- 
cleon isovector magnetic moment is nearly 5. This term is thus ~25 times larger 
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than if it were of "normal" size. A retardation correction contributing to all electric 
multipoles is moderately important, and is the only one of our terms to which the 
nucleon form factors contribute. Finally, higher-order terms, including the seagull 
contribution, are estimated and shown to be tiny. 
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Appendix A 



We wish to evaluate the contributions of Figures (la)-(lc) to the energy of the 
nth hydrogenic S-state. Because this calculation has been set up before|0], we sketch 
that part of the derivation. 

The nuclear energy and momentum scales are much greater than those of an atom. 
Consequently, such large momenta flow through the photon and electron propagators 
in Fig. (1) that only the shortest-range part of the electron wave functions, |0„(O)p, 
contributes to leading order in Za, the product of the nuclear charge, Z, and the 
fine-structure constant, a. Consequently, the momentum in both photon propagators 
is taken to be q (differences in these momenta lead to higher-order terms in Za). It 
is important to enforce the constraints of gauge invariance||15[| on the nuclear part 



of the interaction (the virtual Compton amplitude), and this is most conveniently 
handled using Coulomb gauge, which isolates the nuclear charge density from the 
transverse parts of the currents (Figs, (la) and (lb)) and seagull (Fig. (Ic)). The 
calculation requires a relativistic treatment of the electron (since q 3> me, the electron 
mass), but a nonrelativistic (i.e., leading-order) treatment of the nucleus suffices. 
We expand the nuclear current, J'^, in powers of 1/M, the nucleon mass, and keep 
no powers higher than linear. One immediate consequence of the latter is the lack 
of (nuclear) momentum dependence in the nuclear charge density, p(x) (unlike the 
current density, J(x)), and the nuclear seagull density ||T3|, i?™''"(x, y), which also lacks 



charge components (i.e., B^^ = for yU = or = 0). Our conventions follow Ref. ^8 
and correspond to natural units {h = c = 1). We will incorporate meson-exchange 
currents where required by gauge invariance, although their numerical contribution 
is ultimately small (see main body of paper). We ignore nuclear recoil corrections 
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(~ l/Mt, the total nucleus mass) in the nuclear operators, but maintain them in 
reduced-mass factors in the atomic basis states. 

The energy shift of the nth hydrogenic S-state due to nuclear polarization is most 
conveniently calculated by performing the contour integral over the time component 
of q^{qo) in the loops of Fig. (1), which leads to0] 



pol 



TT 



l0n(O)P 





E 







+ 



+ 



Ami 



2E + UJN 



iq'^ml 



Eq%E + unY - ml 



(2E + cu^)|(iV|p(q)|0)p 
Eq'^liE + unY - rn?^ 

(2g + LUjv) 

Amlq^{q + uj^Y" 



{Al) 



where q^ = c^,E = \Jq^ + ml, uo^ = Ej^ — Eq, the energy of excitation (relative 
to the ground state) of the Nth nuclear state (which by assumption cannot be the 
ground state). In addition 



p(q) 
J(q) 



c/='xe*'i-Xx), 
rf='a;e*'i-"J(x), 



iA2) 

m 

(AA) 



and "_L" signifies contraction with respect to 5*-' — q'^q^ (i.e., J']_ = J^J^{6'^^ — q^q-'))- 
Gauge invariance of this nuclear Compton amplitude (discussed in Appendix B) re- 
stricts B'^ to the "inelastic" part, B^ (see Eq. (B12)). 

Because nuclear momentum, size, and energy scales are given by Q ~ l/i? ~ 100 
MeV, un ~ Q'^/2M ~ 5 MeV, while mg ~ 0.5 MeV, there are two small dimensionless 
expansion parameters: 6 ~ unR ~ 1/20 and P ~ rrieR ~ 1/200. Consequently, 
we work in configuration space, where expansions in R are easiest, using techniques 



introduced in electromagnetically-induced heavy-ion reactions |]T8[. Similar techniques 
were used in Refs. [|], If. We introduce the inelastic transition densities (squared) 



Pi„(x,y) = (0|pt(x)|iV)(iV|p(y)|0), 
JS-(x,y) = (0|J^t(x)|Ar)(iV|J^-(y)|0) 



{A5) 
(AG) 



and Bll{-x,y) has already been used in Eq. (A4). We are not interested in hyper- 
fine structure and, consequently, the spin average over nuclear ground-state (|0)) 
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azimuthal quantum numbers is assumed, while the sum over these quantum numbers 
for the intermediate states (lA^")) is imphcit in Eq. (Al). Consequently, pm and Jf* are 
space scalars. The plane waves in Eqs. (A2 - A4) can then be extracted, appropriate 
spherical averages over q performed (e.g., e**' *^^"^-' — > sm(qz) /(qz), where z = x — y), 
and the g-integrals in Eq. (Al) finally evaluated. This produces the generic result 

AEpoi = -8aV|0n(O)|' / d^x f d'y ^ [pin(x, y)/A,(^) 



+ Ji^(x,y)J^(2) + J2(x,y)zVJ^(z)] 
+ ^Bg{^, y)K{z) + ^Si^(x, y)z'z^K{z) 



{Al) 



where all of our effort will be devoted to obtaining the polarization structure functions: 
In{z), Jn{z), Jn{z), K{z), and K{z). We will develop general forms and then perform 
expansions in {lunz) and {niez). 



We begin with the charge-charge interaction term, 1^, which is the most compli- 
cated to obtain. All other integrals can be obtained from this one: 



Jo 

/ 

Jo 



dq[{2E + u!n)'-^n] sm{qz) 
q^E[{E + UNy-ml] 
°° dq sin(g^) 



ujnZ Jo q^E 
(7W(^) -7W(^)) 



1 - 



{coN + Ey-ml 



Nominally infrared divergent, that part of I'^iz) does not contribute after the integrals 
over x and y are performed in Eq. (A7) (the nucleus is assumed to be virtually 
excited). Ignoring all of these (ultimately vanishing) terms here and in subsequent 
integrals, we find with /3 = rngZ 



Jo 



dq sm{qz) 
q^ 



mi? Jo 



y 



— K r dx r dy r dp' Ko{p') 
Jo Jo 



1 



2ml 



/3 



d(5'{(5-(5'f K,{(5'), 



where Kq{x) is the modified Bessel function of order zero. The function I^^\z) can 
easily be expressed in terms of the Kin{z), the (nth) repeated integrals of Kq{x) [see 
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Eqs. (9.6.25) and (11.2.8) of Ref. ||3^]. The remaining integral J*^^^ is more challenging 
and requires several steps. We define ^ = UN+me and ^' = UN—rrie, assuming C,' > nie 
for the results reported below (this is easily relaxed), and write 



= ^r^sin(g.) 
2mp Jo q 



-1 



+ 



E{E + i) E{E + i') 



/«(e';^)-/«(e;^) 



2mp 



(AlO) 



Using the field-theory trick 

1 



{All) 



EiE2{E^ + E2) TT Jo (A2 + El){\^ + El) ' 
and Eq. (3.737.3) of Ref. the Fourier transforms can be evaluated in the form 

dX (l-e-^V^^) 



i^'\i:z) 



(A2+e) 

7r 



(A2 + ml) 



2me(^ + me) me{i'^ — ml) 
dx X e~^^ 

1 (x^ + /i2)A/a;2 — ij ' 



°° dxe 



1 xya?^ 



1 



{A12) 



after using partial fractions and z'^X'^ = /^^(x^ — 1), and defining /i^ = ^"^ /ml — 1(> 0). 
The first of the remaining integrals is Kii[j3) = it/2 — Jq dj3'Ko{i3') (use x = cosh(t) 
in Eq. (11.2.10) of Ref. [^). The remaining integral (defined to be /(/5; /i) with the 
minus sign) yields to a further trick; it is the solution of the differential equation 



a/32 



+ f,'iiP) = K',ip) 



(A13) 



subject to the boundary condition, J(oo 
solutions exist: 



iiP;fx) 



I 

H J 13 

-m, 



^ exponentially, for which standard 
rf/?'K^(/5')[cos(/i/5') sin(/i/3) - cos(/i/5) sin(/i/?')] 



+ 



^ ,sin(/i/3) log(/i + + 1) + - cos(/i/3)] 
rf/3'i^o(/5')[sin(/i/3) sin(/i/3') + cos(/i/9) cos(/i/3') 



(A14) 



16 



after integration by parts and rearrangement. Combining terms from Eq. (A12) and 
using the last form of Eq. (A9), we obtain /^^^ = + P), where 



m^^e-K) 1 2e (l-cos(/./3)) + Ysin(;./3)log 



2 _ ^2 



- fd(5'Ko{(5') 
Jo 



sin(/x/3) sin(///3') + cos(/x/3) cos(///3') 



(A15) 



Finally, a relatively simple result is obtained 
1 



2meUJNZ 
( 1 _|_ locr ' 



/(e;/3)-7(e';/3) 



TTZ 

~2A 



^nz^ / , / n ^^3 39\ ^, 5, 



(^16) 



where 7 is Euler's constant and 

UN 



2m^ 



2 \ "/2 



I i2 



log 



12 _ j^2 



On log ( ) + 6„ 



rUe 

4 

H ^ 



{An) 



1 + n log 

\ rUe J 



an log ) + hr, 



Useful values are (ai = 1/2, 61 = —5/12), (oi = 7/8, 61 = and (03 = —1,^*3 

j^). This completes treatment of the charge term. 



The current and seagull terms require special treatment because of the transverse 
projection (i.e., the "±"). The appropriate current-current terms have the generic 
form 



jd\fmN\iM)m 



d'x I d'yJli{^,y) I d'qf{q)\6^^ - 



rf^x/ AUi^(x,y) I d'qe'^^-^fiq) 
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+ ^(x,y)v;v^ 



tqz 



/(?) ■ 



(^18) 



The second q-integral (defined to be Jcg) can be seen from Eq. (1) to be 
-^[In{z) — linim^^o -^Af(-2)], wliile tlie first integral (defined to be J°) is given by 

[—j^^-^z [iNiz] — \imin^,^o In{z)]]. Both forms have an infrared divergence (which 
will ultimately disappear because of gauge invariance) and we write mg 2A in 
divergent logarithms, where A is the equivalent small-gr cutoff in the integrals. It is 
easy to perform the derivatives, take the limits, and perform the subtractions. For 
the sake of brevity, we quote only the power-series forms. 



Ami 



J, 



CG 



log(2A/me) + ALi 



log(2A/me) + ALi - 



-AL, 



6 
20 



+ 0{z'). 



Performing the derivatives in Eq. (A18), and matching to Eq. (A7), we find 

1 



J, 



N 



6mj 00 N 



log(2A/me) + AL, - ^AL, 

o 



and 



J, 



N 



The seagull integrals are straightforward variants of Eq. (A9) . We find 
and 



+ 1^(7 + logW2)- 77/60) +0(.'), 



1 
120 



(7 + log(/5/2)-31/30) + O(;.2)_ 



(A19) 
{A20) 

{A21) 
{A22) 

{A23) 
(A24) 



We have kept terms ~ i?^ with the charge densities, since the term of 0(1) leads 
to a vanishing result and therefore the leading order ~ R'^. The currents themselves 
are of leading order R"^, so we have kept terms through i?^ (i.e., of relative order R^), 
as in the charge-density case. 

With hindsight, we can identify the dominant terms in the expansion, which are 
thoroughly discussed in the main text. We note the appearance of log(2;) terms and a 
z^ term. If the Fourier transforms of p and J could be term-by-term expanded in g^, 
these terms would be absent. The integrals diverge at some order, signaling this with 
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"non-analytic" terms in 7? (e.g.,z^ = (z'^)^^'^). Such terms play an important role in 
modern effective field theories[^]. In small terms that we will only estimate, charge 
terms and z"^ seagull terms, we will replace log(z) and log(ti;Ar) by average (constant) 
values, log(^) and log{E), respectively, in order to obtain tractable expressions for 
estimation. 



A tedious application of Cartesian moments (discussed in Appendix B) after ex- 
panding z'^ and z'^ in powers of x and y leads to 

AEn = -8aW\(Pnm^[A + B + C], 



{A25) 



where the charge-charge contribution is 



A 



E 



mm? 



[1 + log(2cujv/me) + --^(log(2cjjv/me) - 5/6)] 



+ ^Jd'xJ rf3y|x-y|3p,,(x,y)+ X;^^(7 + log(^^)-^)[KiV|Q'^ |0)| 



+ -\{N\r'\0)\' - 2{0\B\N) . {N\Om , 
while the current-current contribution is 



(A26) 



B 



+ 



E 

7ml 



( ^ log(2A/m.) + log(2c.^/m.) - 5/6 



12uj 



N 



mt 



N 



20i 



- — ((0|N|iV) • (Ar|D|0) - (0|D|Ar) ■ (Ar|N|0)) 



{A27) 



where Q^^ is the traceless (E2) part of Q'-', while the seagull term is 

5S(x,y) 



C = Jd^xJ 



I2ml 



+ 



+ 



(7 + log(/?/2)-4/3) 



240 
1 

30 



log(2A/me) 

-(g^(x,y)x-y-ij^xV) 
12 



[7 + log(/5/2)] [BliMf - -X ■ y) + 5^i?(-2yy - -xY + 1/V 
(5?*(154y2 - 54x • y) + B^{-mx'y^ - Q2yY + 31y*x-''; 



(A28) 
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where we have separated out a special seagull term (the second term in brackets) that 
has the multipole character of Ml^ and defines the diamagnetic susceptibility[^]: 



Xz>(x, y) = --(5S(x, y) X ■ y - 5?^(x, y) xY] 



(A29) 



We note the small coefficients (~ 1/240) of the last of the current and seagull terms. 
Similar terms that arise from the charge are an order of magnitude larger, as we 
will see through the use of the gauge sum rules of Appendix B. The former terms 
are estimated in the main body of the paper and will prove to be entirely negligible. 
Consequently, we will ignore those current and seagull terms in what follows. We 
also note that the potentially large (infrared) factor of log(me) cancels in these two 
higher-order terms, except for the coefficient of xd- 

Relation (B14) shows that the infrared- divergent terms ~ log(2A/me) cancel iden- 
tically. Using the relations (B15 - B19), we can reexpress the last term in (A26) in 
terms of seagull operators. We finally obtain 



A + B + C 



E 



|(iV|D|0)P 



m 



1 + log i2ujN/me) + —\ (log (2u;Ar/me) - 5/6) 



N 



71 

24 



X 



+ 



J £x J £y |x - y I Vin (x, y) + ^ (7 + log (E 
d'x I d'y (x ■ y - + Bii (-2yy + xY + v'x^ 



39 \ 
20/ 



E 



log(2u;Ar/me) - 5/6 + 



7ml 



log {2ujN/me) - 




^KiVj^,l„g,2.„M)-l/6) 



iv^o 12^^ 



+ 



d'x / d^'y (7 + log {(3/2) - 4/3) xd (x, y) 



(v430) 



It is shown in the main text that only the D, ^^ and z'^ terms contribute at the part 
per thousand level. We have enclosed separately in large brackets the charge-charge, 
current-current, and seagull terms. The last of the charge-charge terms can also be 
rewritten in a form that allows estimation of interaction-dependent (potential) terms: 



AA=^Jd'xJ d'y (7 + log [Ez) - ^) (0| [p(x), [H,, p(y)] ] |0) 



(A31) 
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Appendix B 



In this appendix we perform a Cartesian multipole decomposition of the currents 
and of the virtual nuclear Compton amplitude [|15], [1^. This is unconventional, but 
affords us the easiest mechanism to impose the constraints of gauge invariance. The 
Compton amplitude so decomposed yields gauge sum rules that express the total 
content of gauge invariance in the long- wavelength limit |jl5|, 0. These constraints 
will be imposed on the results of Appendix A. 

The nuclear current is conserved, or 

V-J(x) = -z[ifo,p(x)] , (51) 
where Hq is the internal Hamiltonian (no recoil), and this leads to Siegert's Theorem 



in the long- wavelength limit [16 



/ 



where 



J rf^xxV- J(x) =i[Ho,'D] 
D = y (i^xxp(x) , 



(52) 
(S3) 



thus removing the explicit effect of interaction currents (meson-exchange currents) 
from the electromagnetic-interaction operator in nonrelativistic order (those currents 
are implicitly present in Hq in the last form of Eq. (B2)). We can extend this result 
by expanding Eqs. (A2) and (A3) in a power series in (g*a;*), arranging the Carte- 
sian indices of the x and J according to representations of the permutation group: 
symmetric, antisymmetric, and mixed symmetry. Just as Eq. (B2) allows the zeroth 
moment of J(x) to be equated to the (time derivative of the) first moment of the 
charge density, all symmetric moments of J are determined by that density through 
current conservation. The other-symmetry moments are model dependent and solely 
dependent on the magnetic-moment density, i-i{x.). One finds [p!6| 



p(q) 



d^x p(x) e^'i-^ ^Z + tq-B- -Q'^q'q^ + ■ 



m 



J(q) 



+ 



d^xJM e'"^ '- 



q^O qq O 



30 

g^N qq ■ N 



15 



«q X - 2 



+ 



(55) 
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where 

= J £x^x^p{^) , (56) 



tJi=]^j d\[^x J(x)] , (58) 

N = i y" rf^x [x X (x X J(x))] . (59) 

The first five terms in Eq. (B5) define unretarded-El (D) and retarded-El (O and N) 
interactions, while ^t is the magnetic-dipole operator, and Q*-' the electric quadrupole 
tensor, which generates E2 (via the traceless E2 tensor, Q*"*) and EO (via the trace of 
Q*-') operators. Terms in J proportional to q will vanish in our case because of the use 
of Coulomb gauge. The large contribution of meson currents makes it convenient to 
use this decomposition and eliminate as much model dependence as possible. To the 
order that we are working this model dependence resides in ^ (Ml) and N (retarded 
El), and Eq. (B5) completely and uniquely summarizes the constraints of gauge 
invariance to this order. Moments of the charge and current densities can be obtained 
by taking derivatives with respect to q. 

We can also develop the constraints of gauge invariance for the Compton ampli- 
tude. This is performed in Ref. [0. Replacing q in Eq. (B4) by q2 and q in (B5) by 
qi, the gauge-invariance constraint is 

[J-(qO, p(q2)] = -||p(qi)p(q2) + ^^'^^'"'(qi, q2) • (^10) 

We choose for convenience to divide the nuclear Compton amplitude into two sep- 
arate gauge-invariant parts: elastic and inelastic. The elastic part is the Compton 
amphtude for a point particle of mass, Mt, and charge, Z, multiplied by two factors 
of the nuclear ground-state charge form factor. This requires a seagull operator 

^eKx,y) = ^Po(x)po(y) (511) 

for gauge invariance, while the "inelastic" seagull is then given by 

5S(x,y)=5^^(x,y)-5:](x,y), (512) 
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where po is the nuclear ground-state charge density, normahzed to / d^x po{x) = Z. 
The inelastic amplitude will be gauge invariant if the full amplitude is. With this 
definition, the gauge invariance constraint (BIO) becomes 

[^™(qi),p(q2)] = g2'An"^(qi,q2) - ^ (p(qi)p(q2) - Po(gi)po(g2)) , (S13) 

where a spin-averaged ground-state expectation value is implied. Expanding this in 
powers of qi and q2 leads to gauge sum rules. 

The simplest sum rule results from a single derivative with respect to q^2'. 

[[Ho,-D],--D] = -Jd'xl d'yBi^i^,y). (514) 
One q2 and two qi derivatives produce 



[[^^o,0],-D] 



d'x d'y \E^ + 2B^y^ 



1 



-I [N, ■ D] = - y d'x j dS \Bl y' - Bli y^ 
while one qi and two q2 derivatives generate 

Ho, Q''] , Q''] =-2 [ d'x [ d'y \Bii x ■ y + 5?^ x'y^] + 



D2 



(515) 



(516) 



(517) 



2 d^x d^y 



5^x-y + 5^(a;V-3l/V; 



-A I d'x I d^yBiiy'x^ +4^. 
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+ (518) 

(519) 



In accordance with our earlier discussion, recoil terms such as the last term in Eqs. 
(B17-B19) should be dropped. Other relations are possible, but are not needed. 

The gauge sum rules derived above are rigorous in the nonrelativistic limit. They 
incorporate meson-exchange currents in both currents and seagulls. It is expected 
that such currents could alter the impulse-approximation seagull by a factor of two, 
based on numerical studies|2^ of Eq. (B14) reexpressed in its usual (Thomas-Reiche- 
Kuhn or f-sum rule) form: 



2Y^ujn\{N\B\0)\' = J d'x J d% 5^^(x,y). 



(520) 
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